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ANALYSIS OF MAXIMUM RANGE TRAJECTORIES 
FOR ROCKET-PROPELLED LUNAR FLYING VEHICLES IN 
A UNIFORM GRAVITATIONAL FIELD 

By A. Gary Childs and Ernest S. Armstrong 
Langley Research Center 

SUMMARY 

The problem of obtaining maximal-range soft landing trajectories for a thrust- 
limited rocket having a prescribed amount of fuel under the assumptions of planar flight 
subject to uniform gravitational acceleration and negligible aerodynamic forces is ana- 
lyzed through the Pontryagin maximum principle and the Davidon method of minimization. 
Bounds on the rocket thrust magnitude are treated parametrically. Nondimensional 
curves are given for calculating extremal trajectories for a wide range of vehicle designs. 
Least-square approximate equations linear in the ratio of final mass to initial mass are 
presented for control parameters which may be used to determine extremal control time 
histories completely. These equations could therefore be useful in analytical design 
studies, for example, those for which data are needed on maximum attitude angle, maxi- 
mum attitude rates, and so forth. 

For zero lower thrust bound cases in which the extremal range is short because of 
the use of small amounts of fuel, the solutions approach impulsive results. Also, for 
zero lower thrust bound cases having short burn periods, it is shown that the results of 
the maximum principle for variable thrust direction may be closely approximated by use 
of a constant thrust direction over the thrusting intervals. 

INTRODUCTION 

Rocket flight may become a very practical method for traveling over the lunar sur- 
face. Numerous rocket -propelled flying vehicles have been proposed for transportation 
in the vicinity of the moon. A study of such vehicles is presented in reference 1. In the 
design and use of these vehicles, it is desirable to know the longest range attainable by a 
given rocket for a prescribed amount of fuel. Accordingly, this report analyzes thrust- 
limited rocket trajectories which maximize the range traveled for a given amount of fuel 
and negligible aerodynamic forces. Since men or fragile instruments may be onboard 
the rocket, a soft landing constraint is imposed on the trajectories; that is, the final 
velocity and altitude of the rocket must be close to zero. The study considers planar 


trajectories which have sufficiently low altitude and short range to insure that gravita- 
tional acceleration may be assumed to be constant and the curvature of the moon’s surface 
may be neglected. The dynamic equations are nondimensionalized and the optimal thrust 
control law is determined by use of the Pontryagin maximum principle. Bounds on the 
thrust magnitude are treated parametrically. 

Problem formulations such as described, that is, problems concerning a thrust- 
limited rocket with the simplifying assumptions of planar flight subject to constant gravi- 
tational acceleration and negligible aerodynamic forces, form a fundamental and important 
class in the study of trajectory analysis. Trajectory optimization problems - such as, 
thrust law for minimum fuel, maximum range, or minimum time — using this model have 
been frequently discussed in the literature. (See refs. 2 to 8.) However, the purpose of 
these discussions has, in general, been only to give qualitative results. The purpose of 
this report is to give quantitative results which will be useful in guidance studies. 

The analysis presented may be considered a generalization of the similar quantita- 
tive study by Manci (ref. 9) who considers a minimum fuel-fixed range problem with the 
same dynamic model and the same assumptions as this study but with a zero lower thrust 
bound. The generalization follows from the parametric treatment of thrust bound (which 
allows zero lower thrust bound as a special case) and the duality that exists between the 
minimum fuel-fixed range problem and the fixed-fuel— maximum-range problem. The 
analysis presented also differs from reference 9 in that, for zero lower thrust bound 
cases with short thrusting periods, comparisons of optimal performance are made between 
solutions in which the thrust direction is held fixed over the thrusting intervals and solu- 
tions in which the thrust direction is allowed to vary. 

SYMBOLS 
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effective exhaust speed 
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D,E,F,M,N arbitrary constants in equation (15) 

f = w i[^(tf)] 2 + w 2 [y( _t f)] 2 + w 3 [y(^)] 

f s = Po CO JL" + P 4 - Pc (switching function) 

f (t“) value of f s obtained from lim f s (x) 

x— t 
x<t 

g constant gravitational acceleration 

H = -P 0 x + Pl x + P + P 3* + P 4(- I -f Ll " *) 

m mass 

p = iog e L m U 

\ m o - H/ 

Pontryagin auxiliary variables (i = 0, 1, 

q - ^J^r) 

\m 0 - HJ 

s dummy variable 

T thrust 

t time 

tj_,t 2 switching times 

Wf final weight 

W 0 initial weight 

Wi weights for f (i = 1, . . 3) 

x downrange distance 


y 


altitude 



“i = (Pi+l (0) ) n (i = *’ ■ ■ •’ 3) 



Subscripts: 
f final 


max maximum 

min minimum 

n normalized (divided by -p Q ) 

o initial 

Operations: 

(") ( ) nondimensionalized 

(•) ( ) differentiated with respect to time (that is, with respect to t 

or t, whichever is appropriate) 

(") ( ) twice differentiated with respect to dimensional or nondimensional time 

DYNAMIC EQUATIONS 

The coordinate system chosen for the problem is given in figure 1. Downrange 
distance is x, altitude is y, the thrust magnitude is T (T m i n = T ^ T max ), and the 
thrust angle is 6 . The dynamic equations for the stated problem are based on the 
assumption that the rocket is to be launched from rest at t = 0 and is required to per- 
form a soft landing at t = tj, and are (ref. 10): 


4 


" _ T cos 9 
m 


T sin 6 

y = -s— - g 


\ 

(x(0) = 0; x(0) = 0; x(t f ) = o) 

(y(0) = 0; y(0) = 0; y(tf) = 0; y(t f ) = o) \ (1) 

(m(0) = m Q ; m(t f ) = mf) 


known initial mass 


known terminal mass 


effective exhaust speed of rockets 


g constant gravitational acceleration 

and x(t f ) is to be maximized. As in reference 9, equation (1) may be nondimensional - 
ized by 

Length = Length 


Mass = Mass 


Tmax 


Thrust = 


Thrust 


Time = Time! 


The bars indicate nondimensional quantities. The nondimensional equations corre- 
sponding to equation (1) are 


£ _ T cos e 


f = T 6 - 1 
m 


m = -T 


(x(0) = 0; x(0) = 0; x(tf) = o) 

(y(0) = 0; y(0) = 0; y(tf) = 0; y(tf) = o) > (2) 

(in(0) = in 0 ; m(tf) = inf) 


IL 


where 


Tmin 


T 


max 


ST ^ 1 


The 0) and (") notations are used for differentiation with respect to nondimensional time 
as well as with respect to dimensional time. Only solutions with positive altitude are 
considered. 


OPTIMAL CONTROL 


Equations of Optimal Motion 

Pontryagin’ s maximum principle (ref. 11) was applied to the foregoing maximum 
range problem with the stated equations and boundary conditions. The maximum princi- 
ple is a necessary condition which the control of a problem must satisfy in order to be 
optimal. The controls of the problem are the thrust magnitude T and attitude 6. The 
quantity to be optimized is the range x(tf). The maximum principle states that the con- 
trols which minimize the negative of the range (that is, maximize the range) will be the 
same as those which maximize 


H = - 


p Q x + p x x + p 2 | 


Tcos e Vp 3 yH-P4Ff^-i)- p 5 T 


m 


m 


for every t with the auxiliary variables p i (i = 0, 1, . . ., 5) given by 


p 0 = ° 

i>! = 0 
p 2 = p 0 ” P 1 


( p 0 = °) 
(Pl(tf) = o) 


(3) 

(4) 

(5) 


p 3 = ° 

P 4 = -P 3 


( 6 ) 

(7) 


Ps = 


m 


■^2 cos 6 + p 4 sin dj 


(8) 


The principle also states that when H has been maximized, it will be zero for all 
The part of H that involves T is 


™ ( cos 6 , „ sin 0 \ 

T b— + P 4-ET' P 5) Ws 


6 



where f s is called a "switching function." It is clear that no matter what value 
assumes, the T which maximizes H will be (ref. 11) 



f s (t ) > (T 

f s (t)< o r 


6 


(9) 


For t > 0 so that fs(t) = 0, T(t) takes the value given by equation (9) determined from 
the sign of f s (t ). 

It will be subsequently demonstrated that f s (0) > 0 and singular arcs (where 
f s (t) = 0 over a nonzero time interval) cannot occur. Thus the maximum principle 
asserts that the engine thrust is either fully on (T = T max ) or at its lowest value 
(T = T m i n ). The zeros of f s (t) corresponding to transitions between T = 1 (T - T max ) 
and T = T m in/T’max (T = T m j.n) are called "switching times." 


In order to determine the thrust attitude, that part of 
The terms containing 6 are 


T_ 

m 


^?2 cos 6 + p^ sin 0 'j 


H 


involving 


6 is examined. 


The terms in parentheses may be considered as the inner product 


* 2 \ 


''cos 0\ 
sin 6 , 


To maximize this product it is necessary to choose the unit vector (^9^ to have the 


direction of ( 1; thus, 


P2 


From equations (4) to (7), 



PlCO = Pi(tf) = 0 

P 2 ( t) = p 2 (0) + p Q t 


(10) 
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p 3 (t) = p 3 (°) 

P 4 (t) = P 4 (0) - p 3 (0)t 


From equation (3) Pq is a nonpositive constant. If Pq were zero, then p 2 (t) would 
be constant and cos 6 would have the same sign for all t. However, physically it is 
known that cos 6 must be positive along some portion of the trajectory (in order to 
achieve positive range) and negative along the terminal portion (in order to achieve a 
soft landing). Therefore, p Q must not be zero and p 2 (0) > 0. Since only positive 
altitude solutions are considered 


sin 0(0) > 0 


or 


P 4 (0) 

^p 2 2 (°) + P 4 2 (0) 


whereby p 4 (0) > 0. If p 2 (t), P 3 (t), and p 4 (t) are divided by -p Q , neither e(t) nor 
the switching times are changed and the equations for P 2 (t) and P 4 (t) may be 
expressed in terms of three rather than four constants: 


where 


P2,n(t) =Q! l " 1 

H,rfi)= a 3~ a 2 l 


a l = 


P 2 (0) 

~^0 

p 3 (0) 


a 9 = 

2 -p 


0 

P 4 (0) 


Of Q = 

3 -p 


0 


( 11 ) 

(12) 


An n subscript means that the equations have been normalized by -Pq. 

Next consider the determination of the switching times. Differentiating f s ^ n by 
using equations (2), (8), (10), (11), and (12) gives 
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x s ,n 


4(p 2) n COS 0+P 4 , n Sin 6 ) 

m(t) 


-(p 2) n co s0 + P 4 sine)-^-p n 

m 


d J 2 2 

ifK" +I v 


m(t) 


with 



m (t)J<p(t) 


A = 
B = 

C = 


1 + “ 2 2 

-2(“i + « 2 “3) 

“i 2 + “s 2 


(13) 


and 

cp(t) = At 2 + B t + C 

Since the numerator of f s>n (t) is linear in t, f s n (t) can have at most one zero. By 

Rolle's theorem from elementary calculus, f s n (t) can have at most two zeros. Since 
A > 0, singular arcs cannot be optimal. 

The condition H = 0 at t = 0 gives 


T(0)f Sjn (0) - a 3 = 0 

and implies 

^s,n(0) = “3 (14) 

Thus f S; n(t) is determined by the differential equation (13) with initial condition (14). 

Suppose now that a set (i = 1, 2, and 3) is found which yields a trajectory having 
positive altitude and satisfying the fuel constraint 

m(tf) = m f 
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L 


A set 04 (i = 1, 2, and 3) satisfying the positive altitude condition may be found by 
either trial and error or by the use of impulsive results. The fuel constraint can be 
satisfied by using 

m(t f ) = m f 

to determine tf. The first switching time tj (if it should occur) may be sought by 
use of differential equation (13) for the switching function f Sjn (t). With T = 1 and 
m(t) = m 0 - t, equation (13) can be integrated in closed form over [Ojtjj by making use 
of the indefinite integral (ref. 12) 


Ds 2 + Es + F dg _ _1_ (jo \Jqj(sj + 
(Ms + N) \j<p( s) M C 


2A 


(* -w) 


BD 


2A\/A 



FM 2 - NME + DN 2 


m 2 i/an 2 + CM 2 - BMN 


loge 


/m^CM - BN + (BM - 2AN)s - 2 \J( AN 2 + CM 2 

- BMN)<p(s)J S 

Ms + N 

J 


(M * 0) - 


(15) 


where cp( s) = As* + Bs + C. This operation yields 


ls,n(ll) _ 


y + Am 0 

^ . . log. 




2C + Bm 0 + (B + 2Am 0 )ti + 2^>(m 0 ^^t^ 


m 0 - t! 


2C + Bm Q + 2]Jcp^m 0 jC 


m r 


+ 2Atj + 2 ^Acp^t-^ 

B + 2iAC 

where f s n (0) = has been employed. Then tq is the first positive real root of 

in the interval [o,m 0 - m f j which may be found with Newton's method from 
numerical analysis. If no such root occurs in the interval jo,in 0 - m£j, the rocket con- 
tinuously burns with T = 1 until the allotted fuel is consumed. 

When ti occurs, a second root t£ § ti may exist. For the T m i n = 0 case, 
t2 has a simple relation to ti. Since mass is constant over (ti , t^] , integration of 
equation (13) over (t-^trjj yields 



(16) 
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f 2 ,n CO 

+ P4,n 

Ct) 

m 

(h) 



tl 

and f s,n(tl) = f s,n(t2) = 0 S ives 

P2,n( _t 2) + P4,n( _t 2) = P 2 ,n( _t l) + P4,n( _t l) 

Substituting for p 2 and p 4 n from equations (11) and (12) gives, for t 2 different 
from tj_, 

- B - 2 ( a l + a 2 a 3) - 

A 1 - O 1 


l 2 A 

From to - t]_ > 0 and equation (17), 


h 


1 + a 2 2 


a l + a 2 a 3 
1 + a 2 2 


(17) 


This inequality is useful for applying Newton's method to equation (16). For T min # 0, 
equation (15) may again be applied to integrate equation (14) (with T = T m i n /T max ) with 

m(t)= m(t]) - ^^(t - tj) over (t 1 ,t 2 ] 

' ' 1 max 


to obtain 


£ s,n(‘2) - f s,n(‘l) + 


BM - 2 AN 


2m Van 2 + CM 2 - BMN 


X log € 


2CM 2 - BMN + (BM 2 - 2AMN, 

s 


2M^(AN 2 + CM 2 - BMN’ 

<p(}2) 


m 





2CM 2 - BMN + (BM 2 - 2AMN)t 1 - 

2M^(AN 2 + CM 2 - BMNj 

<p{h) 


ml 




J 


+ ^ log 
M e 


B + 2At 2 + 2\jA(p(t2^ 


B + 2At 4 + 2\jA(p(l.i 'j 


(18) 
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where 


i,n(^l) - 


N 


= ( m o " ti) + 


■ min t 

; *1 

■max 


M = - 


mm 


■ max 


By Newton's method, t 2 can be found as the real root, t 2 > tj, of 

*s,n(t2) = 0 

with f s n(^) gi yen by equation (18). Trajectories in which tj and t 2 exist are com- 
posed of three subarcs - initial and final subarcs with T = 1 and an intermediate sub- 
arc with T ~ T m :m/T max- 

For positive altitude trajectories satisfying the fuel constraint, a single root of 
equation (16) can occur only when f s n (t) has a minimum at t-p In this event the tra- 
jectory is composed of one arc obtained by using T(t) = 1 for tcjo,tf| where 
tf = m 0 - m f . 

For trajectories composed of three arcs, m(t) has the form: 



m(t) = m G - t 

(tc[o,tij) 

(19) 


S(t)=m(t 1 )- (t - tj 

' 1 max s ' 

(tc(ti,t 2 ]) 

(20) 


m(t) = m(t 2 ) - (t - t 2 ) 

(tc(t 2 , 

(21) 

so that 

m f = m(tf) = t 2 + (m Q - t x ) - Tmax (t 2 " h) ~ 

tf 


The final time 

is then a function of the known quantities m 0 , m^, 

ti, and t 2 by the 

equation 

t f = (m 0 - m f ) + (t 2 - tj) - (t 2 - tj) 

1 max ' 


(22) 

With m ( 

j) given by equations (19) to (21) and Pg n (t) and 

p 4 n (t) § iven b y 



equations (11) and (12), the equations of motion for te£o,tfj have solutions 


12 


( 23 ) 


where 


, x , \ . pt (t - s)t/’q ! 1 - s') 

x(t) = x(to) + (t - to)x(t 0 ) + - 1 ds 

t 0 m(s)V«p(s) 

¥rr^ ds 

J to m(s)VMs) 


y(0 = y(t Q ) +(t - t 0 )y(t 0 ) + (t - s) 


T ( a 3~ Q! 2 s ) t 
m(s)^(s) 


ds 


J to m(s)V<p(s) 


ds 


(24) 


(25) 


(26) 


t Q - 0 
T = 1 

m(s) = m Q - 



(27) 


and 





_ T 

T = - 


min 


max 


>tc(ti,t 2 ] 


m(s) = m Q - 


min 


■ max 


( B - ‘i) 


(28) 


^-o _ ^2 


T = 1 


tc(t 2 ,"t f ] 


m(s) = m Q - ti - ~HL(t 2 - tj) - (s - t 2 ) 


max 


(29) 


For a trajectory possessing three arcs, closed-form solutions of equations (23) to (26) can 
be obtained by the successive application of equation (15). For a single maximal thrust 
arc, only equation (15) and equations (23) to (26) with 
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m(t) = m 0 - t 

are needed. 

The Optimal Control Law 

The optimal control law for maneuvering the vehicle described by equation (1) from 
launch to a soft landing while using a given amount of fuel and maximizing the range is 
given by equations (9) and (10). Of course, synthesis of the control law requires equa- 
tions (11) to (29). For positive altitude, the extremal trajectories either are composed 
of one maximal thrust arc or three arcs generated by the thrust sequence: 

(1) A maximum thrust period from t = 0 to t = tf (the first switching time) 

(2) A minimal thrust period from t = tf to t = t 2 (the second switching time) 

(3) A maximal thrust period from t = t2 to the final time tf. 

Once a set of constants (i = 1, 2, and 3) which generates a trajectory having 
positive altitude and satisfying the fuel constraint has been found, a procedure has been 
indicated for constructing the corresponding trajectory. Zeros of equations (16) and (18) 
(obtained numerically) provide switching times. Either equation (22) or tf = m 0 - mf 
determines the final time. Once this information has been given, closed-form solutions 
of equations (23) to (26) which completely specify the trajectory can be obtained with the 
aid of equation (15). There is however no guarantee that the given trajectory will satisfy 
the soft landing conditions. In general, it should only possess positive altitude and satisfy 
the fuel constraint. The next section details a procedure for correcting the given con- 
stants af (i = 1, 2, and 3) to meet the soft landing conditions. 

Determination of Corrected Constants 

The integrated form of equations (24) to (26) evaluated at tf yields equations for 
y(tf), x(tf), and y(tf) explicitly in terms of the a\ (i = 1, 2, and 3). Setting these 
equations equal to zero gives a system of equations in aq which satisfy the soft landing 
conditions. This system can be solved by minimizing to zero with respect to the 
aq (i = 1, 2, and 3) the function 

f (« i,« 2 ’ “ 3 ) = w i *{tf)] + w 2 ypf) + w 3 j(tf) (3°) 

where w^, w 2 , and Wg are positive weights. Constants oq (i = 1, 2, and 3) which 
cause f( a i > a 2’ a 3 ) to vanish and generate trajectories satisfying the fuel constraint are 
the type desired. 
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In the sections to follow, examples of the results obtainable from the foregoing anal- 
ysis are presented. In all cases it was found that an initial set of aq (i = 1, 2, and 3) 
satisfying the altitude and fuel constraint could be easily obtained. The value of 

given by equation (30), with weights w-^ = 10"^, Wg = 10^, and wg = 10“3, 
was iteratively minimized by use of the Langley Research Center Control Data 6600 digi- 
tal computer and the Stewart modification of the Davidon method of minimization. (See 
refs. 13, 14, and 15.) Rapid convergence was obtained to a set of constants cq (i = 1, 2, 
and 3) satisfying all the required conditions. A representative computing time is 20 sec- 
onds. A listing of the program employed is available from the authors at the Langley 
Research Center. 


RESULTS 

The results obtained are of two types: nondimensional and dimensional. Nondimen- 
sional results are completely general with respect to all g, c, and T max> Solutions 
for a range of payload fractions mf/m 0 with initial mass as a parameter are found for 
Tmin = 0 and T m i n - 0.1T max . Although both one arc and three subarc trajectories 
may be extremal, no one arc trajectories were found which satisfied all the necessary 
conditions; only three subarc trajectories are presented. 

Nondimensional Results 

Nondimensional results are given in figures 2 to 8. Least-squares approximations 
are made for the data points of figures 2 to 5. The equations for these approximations 
are presented in appendix A. 

Figure 2(a) shows the nondimensional range achievable with an extremal trajectory 
as a function of payload fraction for three values of nondimensional initial mass (planet 
take-off weight divided by maximum thrust) when T m in = 0. The figure may also be 

used to give the fuel required by a minimum fuel trajectory for a given range and initial 

nff 

mass. The fraction of initial rocket mass used as fuel is 1 - . Note that the range 

m G 

increases with increasing initial thrust-weight ratio. The approximating curves for fig- 
ure 2(a) are quadratic. 

Figure 2(b) is the same as figure 2(a) except that T m i n = 0.1T max . The nondimen- 
sional range is less for given values of m 0 and m^j m Q than when T m i n = 0. Theo- 
retically, this result is implied from the necessary conditions since if T m i n = 0, then 
T = 0 for f s < 0 satisfies the maximum principle whereas T = 0.1T m ax f° r fs < 0 
does not satisfy the maximum principle. Physically, gravity is apparently more efficient 
than rocket thrust for the intermediate period of this type of trajectory. The lessening 
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of nondimensional range is greater for lower m 0 values. This result is reasonable 
since a high thrust-weight ratio (lightweight) vehicle will "waste” a greater percentage 
of its mass during the nonextremal (with respect to the "unconstrained" Tmin = 0 prob- 
lem) intermediate period. The approximating curves for figure 2(b) are quadratic. 

Figure 3(a) shows ti as a function of mf/m 0 for T m i n = 0. If tj is known, 
t2 and tf may be computed from equations (17) and (22). From these results and 
recalling that m G is the inverse of the initial thrust-weight ratio, it is seen that a 
larger thrust-weight ratio will not only produce a longer range but also a shorter (in 
time) trajectory with a longer coast period for the same fuel fraction. The approximating 
curves for figure 3(a) are linear. 

Figure 3(b) gives ti and t 2 as functions of m f /m 0 for T m i n = 0.1T max . 
Comparison with figure 3(a) reveals that for T m i n = 0.1T m ax, the intermediate period 
(*1**2] is longer, but the entire period of flight Jj3,tfJ is shorter. When 
T min = 0-l T max, fuel is burned during the intermediate period and thus the initial and 
final periods are forced to be shorter. It seems reasonable to state that if T m j n is 
close enough to zero, the time lengths of the entire trajectory and of the intermediate 
period would be close to those for T m i n = 0. Thus, if the burn periods are shorter for 
some value of Tmin, and if this value of T m j n is sufficiently near zero, the intermedi- 
ate period should be longer than when T min = 0. Also, for T min sufficiently near 
zero, the entire trajectory should be shorter in time since ail the periods (initial, inter- 
mediate, and final) would be of about the same length as those for T m i n = 0, but fuel 
would be depleted more rapidly since T m i n * 0. The approximating curves for fig- 
ure 3(b) are linear. 

Figures 4(a), 4(b), and 4(c) show, respectively, a^, a^, and <23 as functions of 

mf/m 0 for T m j n = 0. The approximating curves for figure 4 are linear. By using 

an - an t 

these curves and tan 0(t) = — , the approximate attitude may be calculated for 

- t 

any t. The parameters a j, a 2 , and increase with increasing fuel. Since a ^ 

increases more rapidly than does a 3 , the initial thrust angle becomes smaller for 
larger amounts of fuel (and longer ranges). Since a 2 increases, for longer range 
problems 6 will change more rapidly with time along the extremal trajectory. The 
parameters a ^ and a 2 increase with m Q , but a 3 decreases with increasing m 0 
values. Thus, 8 is smaller and changes more rapidly when the rocket vehicle is 
heavier or has a relatively small thrust. 

Figures 5(a), 5(b), and 5(c) give, respectively, a^, a ^ and as functions of 

m f /m 0 with m 0 as a parameter for Tmin - ITmax- A comparison with figures 4(a) 
to 4(c) reveals that when T m i n = 0„lT max , a ^ is greater, a 2 remains about the 
same, and is less. These changes should make cos 6 larger and sin 8 smaller 
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for all values of t. Thus, the values of 6 (or after Q = 90°, the supplement of 6) 
for Tmin = 0.1T max should be smaller than those for T m f n = 0 for all values of t 
for which 8 + 90° once a time translation has been made to make the two 90° points 
coincide. Therefore, the rocket is depending more on the thi’ust and less on gravity for 
obtaining range than in the T m i n = 0 case. This statement is reasonable since the 
rocket can no longer gain as much potential energy from its shorter initial burn period or 
dissipate as much kinetic energy with its shorter final burn period. The approximating 
curves for figure 5 are linear. 

For high thrust-weight ratio (m 0 small) and large values of m-f/ m Q (between 
0.95 and 1) trajectories are produced (for T m f n = 0) which have short ranges and whose 
properties approach those of impulsive cases. Figure 6 (variable attitude) shows non- 
dimensional range as a function of mf/m 0 for m 0 = 0.2. It is a magnification of the 
high mf/m 0 portion of the curve in figure 2(a). Figure 7 shows for m 0 = 0.2 and 
mf/m 0 varying between 0.95 and 1 that maximum altitude increases linearly with range. 
The slope of the line is about 0.225 compared with the analytically determined slope of 
0.250 for the impulsive (m 0 = o) case. (See ref. 9.) Other studies with m Q > 0.2 indi- 
cate that the slope decreases with increasing m 0 . Figure 8 (variable attitude) is a graph 
of tf , t 2 , and tf as functions of mf/m 0 for m Q = 0.2. The coast (t = o) portion is 

seen to be the predominant part of the trajectory. As m 0 ' is further reduced, the 
results approach the impulsive case in which the trajectory consists entirely of a coast 
period with initial and terminal thrust impulses. 

Also, for short-range problems, the rocket attitude changes little during the burn 
periods. Since this condition was noted, a study was performed in which the attitude was 
held constant over each burn period. In the study, the times tf, t 2 , and tf were 
redetermined by taking into account the constant-attitude constraint. Appendix B sum- 
marizes the analysis. 

Figures 6 and 8 (constant attitude), respectively, give nondimensional range and tf, 
t 2 , and tf as functions of mf/m 0 for constant attitude and m 0 = 0.2. These results, 
for mf/m 0 between 0.95 and 0.995, lie so close to those for the variable-attitude case 
that any differences shown on the graphs might be attributable to an error in meeting the 
soft landing conditions in the variable -attitude solutions. The constant-attitude solutions 
satisfy these conditions exactly. Raw data indicate, however, that for constant attitude 
the range is less. No clear trend is obtained for tf , t 2 , and tf from these cases. 

Dimensional Results 

Dimensional results obtained from the computer are presented for two vehicles on 
the lunar surface. The lunar surface gravitational acceleration constant g is approxi- 
mately 5.315 ft/sec2. The two vehicles are: 


17 



(1) A two-man vehicle modeled after a configuration for the lunar-module ascent 
stage with 


T max = 3507.7 lb (15 603.0 N) 
c = 9853.2 ft/sec (3003.3 m/s) 
W Q ~ 10 500 lb (46 706.3 N) 

Wf ~ 8600 lb (38 254.7 N) 


(m Q = 0.49) 

ffiH 


(2) A single-man vehicle proposed in reference 1 with 


Tmax = 275 lb (1223.3 N) 
c = 9660 ft/sec (2944.4 m/s) 

W 0 « 550 lb (2446.5 N) (m 0 = 0.33) 

( EL 
~ = 0 . 
m 0 

Two -m an vehicle.- Figure 9 illustrates dimensional results with T m in = 0. The 
times t-^, t 2 , and tf may be dimensionalized from the values of tj, t 2 , and tf of 
figure 3(a), by using the before-given values of m Q and mf/m 0 , to give the thrust 
history. The first burn period is of longer duration than the second. This difference is 
explained by recognizing that the average mass of the rocket during the second burn 
period is less than it was during the first burn period. 



Figure 9(a) shows the 6 time history for extremal control. For T m f n = 0, 

B is undefined over the coast arc. Points of interest in figure 9(a) may be approximated 
through the nondimensional results. The attitude of the rocket with respect to the near- 
est horizon is 180° - 6 for 6 greater than 90°. The attitude with respect to the near- 
est horizon is seen to be approximately symmetrical for the two burn periods. 


Figures 9(b), 9(c), and 9(d) give the trajectory achieved by the extremal control. 
Any point on the trajectory may be found with equations (23) to (26). The extremal range 

i 

x^ and maximum altitude are about 21 and 34- nautical miles, respectively. The range 
can also be obtained from figure 2. 


One-man vehicle.- Figure 10 illustrates and compares dimensional results for 
Tmin = 0 and T m i n = 0.1T max . The intermediate period is longer, the two burn 
periods are shorter, and the trajectory is shorter (in time) when T m in = O.lTmax 
than when T m f n = 0. 

Figure 10(a) shows the 9 time history for extremal control. For 
Tmin ~ °’lTmax> 6 is g iven by bilinear tangent formula. For T m f n = 0, 6 is 
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undefined over the intermediate period, but for comparative purposes with 

Tmin = 0-lT max , is plotted according to the bilinear formula for the entire period of 

flight. 


Figures 10(b), 10(c), and 10(d) show comparisons for T m i n = 0 and 
Tmin = 0-lT max for the extremal trajectories. The range is shorter and the maximum 
altitude is less when T m i n = O.lTma^. The range for T m i n = 0 is about 23 nautical 


miles whereas for T 


min 


= 0.1T 


max; 


it is about 22 nautical miles. The maximum alti- 


tude for T m j n = 0 is about 4^ nautical miles whereas for T 


min 


0.1T max it is about 


3i nautical miles. For T m i n 


= 0, the horizontal velocity during the intermediate period 
is constant, but when T m i n = 0.1T max , the horizontal velocity climbs to a peak and then 
falls to a value at t 2 comparable with its value at tj. The vehicle reaches a smaller 
maximum and larger minimum vertical velocity during the trajectory when the rocket 
engines are burning during the intermediate period. 


CONCLUDING REMARKS 


The problem of obtaining maximal range soft landing trajectories for a thrust- 
limited rocket having a prescribed amount of fuel under the assumptions of planar flight 
subject to constant gravitational acceleration and negligible aerodynamic forces has been 
analyzed through the Pontryagin maximum principle and the Davidon method of minimiza- 
tion. The problem was solved in nondimensional form with the bounds on thrust magni- 
tude treated parametrically. Nondimensional plots were given to allow calculation of 
extremal trajectories for a wide range of vehicle designs. Least-square approximate 
equations linear in the ratio of final mass to initial mass are presented for control 
parameters which may be used to determine extremal control completely. These equa- 
tions could therefore be useful in analytical design studies, for example, those for which 
data are needed on maximum attitude angle, maximum attitude rates, and so forth. 
Dimensional results were given for lunar flights of a rocket vehicle modeled after the 
ascent stage of the lunar module and a single-man vehicle proposed in NASA CR-365. 

All "best extremal” solutions were found to have similar characteristics. Each 
solution had an initial and terminal maximal thrust period separated by a single period 
at minimal thrust. Thrusting was initially at some positive angle and the vehicle pitched 
up as the initial thrust period continued. After the attitude angle had increased to 90° 
for nonzero lower thrust bounds or after the coast period for zero lower thrust bounds, 
the supplement of the attitude angle descended a path nearly symmetric to that taken pre- 
viously by the attitude angle. For zero lower thrust bound cases in which the extremal 
range was short because of the use of small amounts of fuel, the solutions approached 
impulsive results. Also for zero lower thrust bound cases having short thrusting periods, 
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it was shown that the performance of extremal solutions with variable thrust attitude 
may be closely approximated by solutions in which the thrust attitude is held fixed over 
the burn portions of the trajectory. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., August 11, 1969. 
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APPENDIX A 


LEAST -SQUARES APPROXIMATION 


Least-squares polynomial approximations were made to the data points of fig- 
ures 2 to 5. Quadratic fits were made to the data of figures 2(a) and 2(b) (range as a 
function of mf/m 0 ) since these results showed an obvious quadratic tendency. Linear 
fits were used for the data of figures 3 to 5 since either the results were close to being 
linear or their variations from linearity took on no definite pattern so that a linear fit 
was "as good as any." The data of figures 3 to 5 may be used to determine extremal 
control completely so that the "be'st" linear fits might be useful in the analytic design of 
a simple suboptimal system. The equations for the approximating polynomials are: 


For figure 2(a): 


mf 

Range = 0.30256 - 0.61506 + 

m Q 


mf 

Range = 0.27241 - 0.55611 ~ + 


Range = 0.24394 - 0.50028 — + 

m G 


0.31284 


0.28413 


0.25688 



For figure 2(b): 


mf 

Range = 0.26854 - 0.54531 + 

m Q 


/- \ 2 


/ m f \ 

/ \ 

0.27705 — 

(m 0 = 0.2j 

W 



For figure 3(a): 


mf 

Range = 0.26440 - 0.53946 ~ + 0.27548 

m Q 


\2 



(m 0 = 0.4) 


mf 


Range = 0.24195 - 0.49606 + 0.25463 


m r 


mf 

t X = 0.10634 - 0.10666 ^ 

m Q 

iiif 

t x = 0.21646 - 0.21739 

m Q 


/— \2 


/ mj \ 

(m 0 = 0 . 6 ) 

w 

#1 

11 

0 


< 5 ^ 

0 

11 
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APPENDIX A 


For figure 3(b): 


For figure 4(a): 


For figure 4(b): 


h 


= 0.33599 - 0.33820 


Hlf 

m 0 


tj = 0.07348 - 0.07358 
tj = 0.19178 - 0.19242 
tj = 0.32021 - 0.32205 
t 2 = 0.72527 - 0.72961 
t 2 = 0.67959 - 0.68580 
t 2 = 0.62877 - 0.63763 


IXlf 

m 0 

mjp 

m 0 

m 0 

mf 

m Q 

m Q 


m-f 

m 0 


= 0.41055 - 0.41279 
a 1 = 0.43161 - 0.43463 
a 1 = 0.45551 - 0.45965 


ni£ 

m 0 


mf 

“o 

mf 

xn 0 


mf 

oit\ = 0.04901 - 0.04994 
z m 0 


an = 0.11675 - 0.11761 ~ 
z m o 

mf 

a 0 = 0.25434 - 0.25323 ~ 
z m 0 


(m Q = 

0.6) 

(m Q = 

0.2) 

(m 0 = 

0.4) 

(m 0 = 

0.6) 

(“0 = 

0.2) 

(m 0 = 

0.4) 

(m Q = 

0.6) 


(m Q = 0.2) 
(m 0 = 0.4) 
(m 0 = 0.6) 


(m 0 = 0.2) 
(m 0 = 0.4) 
(m 0 = 0.6) 
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For figure 4(c): 

cn g = 0.38492 
ex 3 = 0.37924 
oi 3 = 0.37934 

For figure 5(a): 

a 1 = 0.39958 
Q! 1 = 0.42820 
= 0.45401 

For figure 5(b): 

« 2 = 0.08132 
a 2 = 0.11943 
a 2 = 0.25498 

For figure 5(c): 

a 3 = 0.29823 
Q!g = 0.35315 
a 3 = 0.37116 


mf 

0.38721 =J- 
m 0 

^m Q = 0.2) 

mf 

0.38253 =J- 
m Q 

11 

0 

mV 

0.38444 =J- 
m Q 

(m 0 = 0.6) 

mf 

0.40178 — 
m o 

11 

0 

mf 

0.43120 =J~ 

m Q 

? 

11 

0 

mV 

0.45809 — 
m o 

'^ 0 ' 

0 

11 

if. 

mV 

0.08458 — 
m G 

(m 0 = 0.2) 

m f 

0.12043 — 
m Q 

O 

II 

if 

mV 

0.25390 zJ- 
m 0 

( m o = 0.6) 

nL 

0.29957 ~ 
m G 

''ST' 

o’ 

11 

jf 

m f 

0.35593 — 
m Q 

0 

II 

if 

in f 

0.37594 — 
m Q 

(m 0 = O.e) 
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APPENDIX B 


CONSTANT -ATTITUDE STUDY 

This section outlines an analysis of the optimal trajectory problem of this report 
with the added assumptions that the optimal trajectory is composed of two maximal thrust 
arcs separated by one coast arc and that the thrust attitude 6 is held constant over the 
thrusting periods. By defining and 0 2 to be > respectively, the fixed attitudes over 
the initial and terminal burn periods, equations (2) may be integrated and the conditions 

x(t f ) = 0 

y( t f ) = 0 

y(tf) = 0 

m(t f ) = m f 


applied to determine the unknowns t 2 , tf, and 0 2 in terms of the single unknown 
tf. This procedure yields 


tf = - 


B ± 




4AC 


2A 


(Bl) 


where 


h - tf - ( m o - m f) + 


sin 0\ 


sin 0 2 = 


P 2 - Q 2 - t f 2 
2Qt f 

P 2 - Q 2 + t f 2 
2Ptf 


m r 


P = log e — 


Q = log p — 


m G - tf 
inf 


m Q - tf 


(B2) 


(B3) 


(B4) 
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and 

A = Q[ti - P(m Q + m f )J + P(m f - m 0 + tf) 

B = PQ(P 2 - Q 2 ) 

C = P 3 |Q(m£ - m G ) + m 0 - m f - tfj + P 2 Qti + p[q 3 (m 0 - m f ) + Q 2 ( mf - m Q + ti)J - Q 3 tj 

The value of tj which maximizes x(tf), if equations (Bl) to (B4) are known, can then 
be found by scanning the allowable values of tf (0 < tf < m 0 - mf) and observing x(tf). 
After determining tj, the complete trajectory may be constructed by use of equa- 
tions (Bl) to (B4) and the solutions of equations (2). 
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Figure 1.- Coordinate system for "flat" surface maximum range problem with thrust components, magnitude T and attitude 9. 
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°80 62 84 86 



m f/ m 0 

T min = °- 1T max- 

Figure 3.- First switching time tp second switching time t 2 , and final time tj as functions of the ratio of final mass to initial 
mass mf/m 0 with nondimensional initial mass Tn 0 as a parameter. 













0 Variable attitude 
□ Constant attitude 


950 .955 .960 .965 .970 .975 .960 .985 .990 .9 

m^/frio 


Figure 8.- First switching time t 1( second switching time t 2 , and final time tf as functions of ratio of final mass to initial mass mf/m 
for very short range trajectories using variable and constant attitude. 7n 0 = 0.2; T m j n = 0. 


CO 





6XI0 4 



I ! I I I I I I 1 J 

0 .5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5XI0 4 


x, m 

(b) Variation of altitude y with downrange distance x. 
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(c) Variation of altitude y with time t. 
Figure 9.- Continued. 
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Time, t, sec 

(d) Variation of horizontal velocity x and vertical velocity y with time t. 
Figure 9.- Concluded. 




Time, t, sec 


(a) Variation of thrust attitude 8 with time t. 

Fiqure 10- Optimal control and trajectories. W 0 ~ 550 lb (2446.5 N); Wf = 450 lb (2001.7 Nl ; T m ax 

c = 9660 ft/sec (2944.4 m/s). 


= 275 lb (1223.3 N); 







x, y, ft/sec 



Time, t, sec 


y, m/sec 


(d) Variation of horizontal velocity x and vertical velocity y with time t. 
Figure 10.- Concluded. 
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